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By considering the local vector current between nucleon states and imposing charge conservation we determine, 
for 0{a) improved Wilson fermions, its renormalisation constant and quark mass improvement coefficient. The 
computation is performed for both quenched and two flavour unquenched fermions. 



1. INTRODUCTION 

Due to the presence of the 'Wilson term' in the 
lattice fermion action for Wilson fermions the dis- 
cretisation errors are 0{a). As the gluon part of 
the action (sum of plaquettes) has only 0{a^) er- 
rors, it is also desirable to achieve this for the 
fermion action. The Symanzik programme^ al- 
lows a systematic reduction of errors to O(a^) 
upon including additional higher dimensional op- 
erators. The (on-shell) action is improved with a 
suitably tuned 'clover' term. However it is also 
necessary to improve each operator separately. 
Much work has been devoted to this topic in 
recent years; here we shall just concentrate on 



the local vector current: V, 



(I) 



In this 

(9) 



case just two additional operators arrig V;i''> and 

^io-dxT^^'^^ are required giving the 0(a) improved 
and renormalised vector current as 



Vi"^'' = Zv{l + am,&y)(U(«) + IzacydMl)) 



with T^^^ = qcTf^^q. The second improvement 
operator only has an effect in non-forward ma- 
trix elements and will not be considered further 
here. The renormalisation constant Zy and im- 
provement coefficient by are functions of the cou- 
pling constant 50 and perturbatively we have, 
[2], to one loop (independently of the presence 
of fermions), 

Zvigo) = 1 - 0.12943.g2 -)-... , 

^'1^(50) = l + 0.15323.g2 + ..., 

but in presently accessible regions of P = G/qq 
there may be considerable deviations. 

2. THE CONSERVED CURRENT 

There is an exact symmetry of the action q 
e~^"q, q e*"g giving via the Noether theorem 
the Ward identity (WI) 



-q 



_dn 

^ dq 



*Talk given by R. Horsley at Lat02, Boston, U.S.A. 
^For a recent review see, for example, [1]. 



where fl is an arbitrary operator, is the lattice 
backward derivative, and jjf"^ the exactly con- 
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served vector current^ {CVC), 

Roughly speaking the RHS of this equation 
counts the number of q and q in operator fl. For 
our numerical results we take O — > BB where B is 
the (unpolarised) nuclcon operator to give, upon 
solving the WI equation, 



{Bit)jt-''\T)Bm 
{B{t)Bm 



<T <t 
t<T <Nt 



where c,-" arc constants 



d) 



(u- 



d) 



and with jump, AR{J\ ' ) = ' — c'2 
1, ie charge conservation. The numerical advan- 
tage of considering u — d\s that the hard to com- 
pute quark line disconnected terms cancel. (For 
the CVC this term vanishes though.) In Fig. 1 
we show this ratio for the conserved vector cur- 
rent. A very good signal is observed (indeed the 
result should be true to machine accuracy). 

3. THE LOCAL CURRENT 

The local vector current {LVC) is not con- 
served on the lattice and so we do not expect 
the jump to bo equal to one. This is shown in the 
RH picture in Fig. 1. We now define the renor- 
malisation constants {Zy, by) by demanding that 
the rcnormaliscd local current has the same be- 
haviour as the conserved current, so that 



Zv{l + amqbv) 



Ai?(yi"-'^): 



Thus upon plotting the data, the intercept gives 
Zv and the gradient .^yfey. {aniq = ^{I/k—I/kc) 
and Kcido) is estimated from roniq oc (roTOpg)^.) 
In Fig. 2 we show quenched results from which 
the intercept and gradient can be found. Other 
alternative non-perturbative determinations have 



^O(a) improvement for non-forward matrix elements re- 
quires, as for the local vector current, an additional oper- 

A)]- 



ator ^ia,^ 



.Tl£(x)+d^Tl£(x- 
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Figure 1. R{.]'f and R{vj" plotted against the 
operator position t for the quenched (uj = 0) data set 
P = 6.0, K = 0.1342 onaMxNx = 16^ x 32 lattice with 



t = 13. Typical fit intervals for 
pairs of vertically dashed lines. 



(u-d) 



are given by the 



been given by the ALPHA collaboration, [3], us- 
ing the Schrodinger functional, the LANL collab- 
oration, [4] using other Ward identities and Mar- 
tinelli et al., [5] by 'mimicking' perturbation the- 
ory. 

As well as quenched data sets (n/ = 0), in col- 
laboration with UKQCD, we have also generated 
unquenched data sets. In this study we use the 
configurations with parameters given in Table 1. 
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Volume 


Trajs. 


Group 


5.20 


0.1342 


le^* 


X 32 


5000 


QCDSF 


5.20 


0.1350 


16-^ 


X 32 


8000 


UKQCD 


5.20 


0.1355 


16^ 


X 32 


8000 


UKQCD 


5.25 


0.1346 


16^ 


X 32 


2000 


QCDSF 


5.25 


0.1352 


16^ 


X 32 


8000 


UKQCD 


5.25 


0.13575 


24^ 


X 48 


1000 


QCDSF 


5.29 


0.1340 


16^ 


X 32 


4000 


UKQCD 


5.29 


0.1350 


16^ 


X 32 


5000 


QCDSF 


5.29 


0.1355 


24^ 


X 48 


2000 


QCDSF 



Table 1 

Data sets used in the unquenched, nj 



■ 2, simulations. 



We now present our results. In Fig. 3 we show 
Zv and in Fig. 4, by- For quenched fermions good 
agreement with other methods is found. 
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Fi gure 2. AR(V'4 ) /or quenched configurations for 
(3 = 6.4, 6.2 and 6.0 (upper set of curves, top to bot- 
tom respectively) and for the unquenched configurations 
for fi = 5.29, 5.25 and 5.20 (lower set of curves). 
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Figure 3. Zy (LVC, filled squares) determined in 
this work for quenched and unquenched 0(a) improved 
fermions. For the quenched case a comparison is made 
with ALPHA, [3], LANL, [4] and RI-MOM, [5]. Pade 
fits are also given for our and the ALPHA results. 

4. CONCLUSIONS 

The method described here reproduces the 
results of other approaches for 0{a) improved 
qTicnchcd fermions. (But one needs to remem- 
ber that Zv definitions can vary by 0{a^) while 
by definitions may vary by 0{a).) For 0{a) im- 
proved unquenched fermions Zy is smaller and by 
larger than for quenched fermions at the same lat- 
tice spacing (roughly a„j.=2(5.25) ~ anj-=o(6.0)). 
Further details and final results will appear in [6]. 
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Figure 4. by (LVC, filled squares) determined in this 
work again for both quenched and unquenched 0(a) im- 
proved fermions. Also shown is the Pade fit from ALPHA, 
[3], and the LANL, [4J results for quenched fermions. 
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